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The quantization of the Friedmann-Robertson- Walker spacetime in the presence of a negative 
cosmological constant was used in a recent paper to conclude that there are solutions that avoid 
singularities (big bang-big crunch) at the quantum level. We show that a proper study of their 
model does not indicate the it prevents the occurrence of singularities at the quantum level, in fact 
the quantum probability of such event is larger than the classical one. Our numerical simulations 
based on the powerful variational sine collocation method (VSCM) also show that the precision of 
the results of that paper is much lower than the 20 significant digits reported by the authors. 

PACS numbers: 98.62.Sb, 04.40.-b, 04.70.Bw 

We review the results presented in a recent paper on Quantum Cosmology by Monerat et al Q], concerning the 
quantization of the Friedmann-Robertson- Walker spacetime. The quantum dynamics of this model is described by 
the Wheeler-De Witt equation, which is obtained by applying the Dirac formalism to the quantization of constrained 
systems. Such equation reads 

T2c& _3fca2+Aa4 ) Vl/(a ' T) = ^c^ Vl/(a ' r) ' (1) 
and it is formally equivalent to the Schrodinger equation for a quartic anharmonic oscillator with a potential given 

by Q 

V e (a) = 36fca 2 - 12Aa 4 . (2) 
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Here a is the scale factor in the FRW metric, A is the (negative) cosmological constant and k denotes the constant 
curvature of the spatial sections. 

The corresponding time-independent Schrodinger equation then reads 

-^+V e (a))r,(a) = 12Er ] {a), (3) 



da 2 

in the notation of Monerat et al [lj]. In what follows we write s = 12 E. 

Our critique of the arguments of Monerat et al [l| is based on both technical and conceptual grounds. The former 
stems from the numerical results obtained in their paper which are below the claimed accuracy. The latter concerns 
the correct interpretation of the model which does not support the conclusion that singularities are avoided at a 
quantum level l|. 

Regarding the technical part. Monerat ct al [lj] obtained their results by means of a technique developed by 



Chhajlany and Malnev |3j, |4[. The application of this method requires the introduction of an extra term in the 
Schrodinger equation, in the form of a scxtic potential with a coefficient obtained from a cubic algebraic equation. 
Due to this fact the energies and wave functions calculated by Monerat et al are only approximations with roughly 
1% error. However, the results reported in Table I of Ref. [lj], which are obtained using the modified potentials given 
in their equations (26), (30) and (32) are claimed to have 20 digits accuracy. Long time ago Fernandez Q noticed 



ti P, u 

terms to reach that precision. 

In order to verify the accuracy of the results we have used the powerful "Variational Sine Collocation Method" 
(VSCM) [a, S|, which allows one to obtain very precise numerical results, both for the energies and wavefunctions, 
and has proved to yield errors that decay exponentially with the number of elements (sine functions) used. 

In what follows we briefly sketch the VSCM. The sine functions are given by 

a (h \ - sin (ttQe - kh)/h) 

SkM = n(x-kh)/h ' (4) 

where k = 0, ±1, ... is an integer and h is the spacing between two contiguous sine functions. The sine functions are 
orthogonal: 

r + oc 



S k (h, x) Si (h, x) dx = h 8 kl . (5) 
Using a collocation scheme, one can express the matrix element of the Hamiltonian operator in the sine "basis" as 



where 



g> + S kl V(kh)} , (6) 



(2) = J -Jh? if k ~ l m 



2 (-I)*-' , f u,j 

are the coefficients obtained from the discretization of the second derivative. Diagonalization of the N x N matrix 
yields a set of N eigenvalues (energies) and eigenvectors (wavefunctions). Notice however, that the diagonalization of 
H requires the specification of the otherwise arbitrary grid spacing h. For a given number TV of sine functions there 
exists an optimal value of h which provides the smallest errors. As shown by Amore and collaborators such 
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TABLE I: Lowest energy levels for the effective potentials with k = 0, ±1. 



Level 


k = 1 


k = 


k = -1 


Ej 


1.5103016760578707333 


0.3404027974289474964 


-10.553138375809812184 


E 2 


3.5509871014722926520 


1.0496018816142096406 


-8.8509608822646650503 


E B 


5.6230931685079239340 


1.9248782225982734499 


-7.2247835560429854525 


Ei 


7.7256531719676264015 


2.9232145294809850623 


-5.6802949920840011967 


E 5 


9.8577817763092210058 


4.0227383298745516995 


-4.2251989178977204001 


E 6 


12.018664366203979043 


5.2097960735587976352 


-2.8707515711793779600 


E r 


14.207548236413247839 


6.4748734068112171253 


-1.6348253113363418720 


E 8 


16.423735233636153025 


7.8108774875618647285 


-0.5364307199618285749 


E 9 


18.666575557925837760 


9.2122719234208330574 


0.481759254997070592451 


Eio 


20.935462499941518896 


10.674588220211030332 


1.583190805382159479497 


En 


23.229827940734217332 


12.194126732921328651 


2.833031712002511140698 


El2 


25.549138478317595566 


13.767762283584094708 


4.213983702088123106309 


Er, 


27.892892073596735435 


15.392811775837602761 


5.707834132050646278314 


En 


30.260615129836301427 


17.066940597587263816 


7.302868225546844668469 


E 15 


32.651859936513780000 


18.788094386112582846 


8.990905434198164180404 


El6 


35.066202421382313215 


20.554447991460572182 


10.76576537713456435413 


E\7 


37.503240164769059826 


22.364366463785281446 


12.62252647249543208197 


Eis 


39.962590638221848703 


24.216374668801875820 


14.55712323817216360861 


-El 9 


42.443889636078437368 


26.109133235245860851 


16.56610649829916492072 


E20 


44.946789873733933173 


28.041419241199542663 


18.64649041909474222481 


E21 


47.470959730597782399 


30.012110508653406479 


20.79564954975016308111 


E22 


50.016082119171135490 


32.020172687758339925 


23.01124646297271199397 


E23 


52.581853464498535137 


34.064648527576543759 


25.29117906504417774373 



optimal value of h can be found by application of the Principle of Minimal Sensitivity (PMS) |8|] to the trace of the 
N x N hamiltonian matrix. Notice that once the hamiltonian matrix has been diagonalized one can also solve the 
time-dependent Schrodinger equation, in terms of the stationary states previously calculated. 

Our table U shows the odd energies for the potentials given in equations (26), (30) and (32) of Ref. [lj using the 
VSCM. The results in this table are correct up to the 20 th digit, and show that the results reported by Monerat 
ct al in their Table I [l| are considerably less precise than they appear to suggest. Just to mention one example, 
their energy E\q corresponding to k = 1 has just 5 correct digits. Although this lack of precision does not affect the 
authors' analysis [lj|, which is in any case bounded by the more severe 1% error due to the use of an effective potential 
rather than the exact one, we believe that it is important to remark this point. 

We wish to stress that the numerical method that we have used is completely general and can be applied to arbitrary 
potentials; for this reason we have decided to perform the numerical simulation and time evolution of the model with 
the exact potential, rather than following the strategy of Monerat et al [lj], which was dictated by their choice of 
numerical method [3, [J. Table ITTl shows the eigenvalues corresponding to the exact potential Eq. (|2J). Particularly in 
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TABLE II: Lowest energy levels for the exact potential with k — 0, ±1. 



Level 


k = 1 


k = 


k = -1 


Ej 


1.5102625384705021518 


0.3364795921008755186 


-21. 79569603947057798056 


E 2 


3.5506472908620463050 


1.0311990503013958552 


-20.39848968823294656013 


E 3 


5.6218937064348010503 


1.8807615814216099958 


-19.01885125967730517826 


E 4 


7.7227778139934265020 


2.8424874930254863943 


-17.65745105165687339198 


E 5 


9.8521942203839334149 


3.8947462112216761884 


- 16.31503013488920386315 


E 6 


12.009138568080285499 


5.0240915609641800702 


-14.99241327017434870429 


E r 


14. 192693387434484198 


6.2211924298081547999 


-13.69052532294706014390 


E 8 


16.402016547726087957 


7.4791208557098708821 


-12.41041248037188872561 


E 9 


18.636331727173164959 


8.7924907241645713076 


-11.15327022078174141020 


Eio 


20.894920472196314023 


10.156971843165374589 


-9.920481041305256697868 


En 


23.177115522582541829 


11.568992820956286787 


-8.713666752230586250096 


El2 


25.482295155850129982 


13.025548149708503746 


-7.534763378250629251892 


Er, 


27.809878360230719403 


14.524066967082505523 


-6.386132854421 764354442 


En 


30.159320687349375640 


16.062320376957969417 


-5.270738293115111860300 


E 15 


32.530110667020528246 


17.638353963148189178 


-4. 192437789021558602284 


Eie 


34.921766690427654626 


19.250437372785090866 


-3.156519423136907030485 


En 


37.333834286304756785 


20.897025823964631886 


-2.170719224552462670164 


Eis 


39.765883728998739663 


22.576730162849598220 


-1.246288323062400948910 


Eig 


42.217507928478425449 


24.288293189747778350 


-0.389996330189812857692 


E20 


44.688320561206965683 


26.030570672865707451 


0.433719867112653240009 


E'21 


47.177954407853964266 


27.802515928224750562 


1.300741393435948191383 


E22 


49.686059869496678687 


29.603167154324565317 


2.243855389930394850590 


E23 


52.212303638550624930 


31.431636924008847119 


3.256001768433707607172 



the case k = — 1, corresponding to a double well we have verified that the results are quite different from the ones 
obtained with the effective potential (compare the fourth columns of tables U and [TTJ) . 

By means of the VSCM we have obtained accurate energies E n and wavefunctions 4> n {a) up to n = 23. Using these 
states we can build a wave packet 

0(o,t) =]TB„0„(a) e—" T (8) 

n 

where the coefficients B n arc taken to be equal to unity as in Ref. Q] and the energies e„ = l2E n . The expected value 
of the scale factor is 

_ f °°a \e(a,T)\2da _ 
W f °°\e(a,T)\*da EjBn? ' W 

where K nn ' = /q°° a 4>n( a ) < t>n'(Q-)da. This expectation value oscillates in time because of the time-dependent terms 
e l (e n -e n ,)t_ Likgwigg the energy of the wave packet is £ wp = ^2 n \B n \ 2 e n / J2 n \B n \ 2 - 



5 



We notice that the wave packet of eq. ([8]) is not yet fully specified, since the eigenfunctions (b n (a) are known up to 
a phase factor which can be chosen arbitrarily. By looking at the coefficients in Table II of [1| and of Tables II- VII 
of we see that these wave functions have been chosen by Monerat and collaborators to be positive in a small 
neighborhood of a = + . The wave packet obtained with this prescription for k = 1 corresponds at t — to a state 
sharply localized around a = 0.1, which quickly spreads at later times according to the Heisenberg principle. The 
reader can look at the left plot in Fig. [TJ where the wave packet at t = and t = 10 is plotted . The wave packet for 
k = — 1 is plotted in the right plot of Fig. [T] 

Fig.l of [lj displays the expectation value of the scale factor for A = —0.1 and k = 1; in this plot it is difficult to 
appreciate the localization of the wave packet, since the curve is compressed by the large time scale considered. In 
figElof the present paper we show the behavior of the expectation value (a) for A = —0.1 and k = ±1 using the same 
choice of phases of Monerat et al. Our result for k = 1 is similar to the one observed by Monerat et al.: the expectation 
value of the scale factor starts at a « 0.1 at t = and then strongly oscillates at later times, in accordance with the 
Heisenberg principle; for k = —1, however, our curve is different from the corresponding curve of Fig. 3 of and (a) 
oscillates around a larger value. The reason for this discrepancy is the poor accuracy of the results of Monerat et al. 
especially for the case k = — 1: the fact that (a) oscillates around a ~ 3 in our figure reflects both the incorrectness 
of the eigenfunctions calculated in [l| and the presence of a minimum of the potential at a = y^f^" = \/l5 (being the 
energy of the wave packet negative). 

The dashed horizontal lines in our figure correspond to the classical expectation value (a) c ; for a particle moving 
in the same potential and with the same energy of the wave packet. The classical probability of finding a given value 
for the scale factor is given by 

P cl {a) = ■ Ncl = (10) 
y/V(a + ) - V(a) 

where N d = l/ J Q a+ ^ y(Q d °_ v(a) and a + is the right inversion point (remember that in our problem there is an infinite 
wall at the origin) given by £ wp = V(a + ). 

Based on the behavior of (a) in their figures 1-3 the authors of Q] state that they "may be confident about all values 
assumed by (a), in particular, that it does not vanish", thus concluding that this proves the absence of a big-crunch 
in the model at the quantum level or using their own words "since the expectation values of the scale factors never 
vanish, we have an initial indication that these models may not have singularities at the quantum level" . 

In our opinion these observations are not correct: in the first place, one concludes that (a) without any 
calculation because it is a consequence of the presence of an infinite wall located at a = and is therefore completely 
independent of the form of the potential. In the second place, and more importantly, the condition (a) does 
not imply the absence of a big crunch and actually holds even at the classical level. As a matter of fact, simple 
mathematical considerations, based on the form of P c j(a), are sufficient to conclude that (a) c / > x . As we have said 
the horizontal curves in Fig. [2] correspond to the classical expectation values which in both cases fall in the region 
of oscillation of the quantum expectation value. For this reason we believe that hardly any conclusion can be drawn 
from the study of (a) performed in Fig. 1-3 of [l|. 



1 The quantum and classical expectation values are obtained as integrals of a positive definite function and therefore it can never vanish. 
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FIG. 1: Wave packets for A = —0.1 and k = 1 (left plot) and k — — 1 (right plot) at two different times, (color online) 




X X 



FIG. 2: Behavior of the expectation value (a) for A = —0.1 and k = 1 (left plot) and k = — 1 (right plot), (color online) 



We also stress that, although the expectation value of the scale factor can never vanish, one can build wave packets 
which at a given time (say t = 0) are localized around an arbitrarily small values, ao- Such wave packets are obtained 
using B n = <j>n{ao) for the first N states of the potential using the completeness of the basis {</>„(a)}: 

N 
n=0 

On the other hand, a more meaningful comparison between the classical and quantum cases can be made by 
calculating the probability of finding a in a region between and a max . Fig. [3] shows the quantum probability that 
the scale factor is smaller than a given value a at three different times, for the wave packet with B n = 1, and k = ±1 
(left and right plot respectively) and also its classical analogue (solid line), using the same phase convention of [l|]. 
We have used a logarithmic scale in the plots to better appreciate the region around a = 0. 

For k = 1 we see that at t = the quantum probability differs strongly from the classical one, since the scale factor 
is localized around a w 0.1, while at later times when the packet spreads, the quantum probability gets closer to its 
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classical 
1 = 
t= 10 
t = 100 
t= 1000 



FIG. 3: Probability of finding a scale factor < a for A = —0.1 and k = 1 (left plot) and k = — 1 (right plot), (color online) 



classical counterpart. At t = only for quite small values of a the classical probability is larger than the quantum 
one, as a consequence of the Dirichlet boundary conditions on the wave functions. 

More dramatic differences are observed for k = — 1: in this case the region around a = is classically forbidden 
and only the quantum probability is nonvanishing around a = 0. 

Concluding, we feel that the question posed in [lj], whether in the quantum model the singularity is avoided, has 
not been addressed properly in that paper : the observation that (a) > is trivial and holds both at classical 
and quantum level; The expectation values calculated for the wave packets considered in [lj oscillate around the 
corresponding classical expectation values and can never vanish neither in classical or in quantum mechanics. The 
strong oscillations observed for the case k = 1 are understood as a consequence of the Heisenberg principle. 

On the other hand, we feel that a better tool to address the question posed by Monerat et al. is to calculate the 
probability of finding a scale factor below a given value (such tool was not considered in [lj ) . Our result show that only 
for very small a the classical probability is larger than the quantum one, although this effect is a mere consequence of 
the Dirichlet boundary conditions, and therefore holds regardless of the binding potential used (in this case a quartic 
potential). 
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